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The one dimensional anisotropic Kondo-necklace model has been studied by several methods. 
It is shown that a mean field approach fails to gain the correct phase diagram for the Ising type 
anisotropy. We then applied the spin wave theory which is justified for the anisotropic case. We have 
derived the phase diagram between the antiferromagnetic long range order and the Kondo singlet 
phases. We have found that the exchange interaction (J) between the itinerant spins and local 
ones enhances the quantum fluctuations around the classical long range antiferromagnetic order 
and finally destroy the ordered phase at the critical value, J c . Moreover, our results show that the 
onset of anisotropy in the XY term of the itinerant interactions develops the antiferromagnetic order 
for J < Jc. This is in agreement with the qualitative feature which we expect from the symmetry 
of the anisotropic XY interaction. We have justified our results by the numerical Lanczos method 
where the structure factor at the antiferromagnetic wave vector diverges as the size of system goes 
to infinity. 
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I. INTRODUCTION 

Quantum phase transitions between Kondo singlet and 
antiferromagnetically ordered states such as found in 
heavy fermion compounds have been attracted much re- 
search interest recently-i^. It is therefore of great impor- 
tance to understand the approach to the quantum critical 
region within suitable theoretical models. The most im- 
portant among them is the Kondo lattice model consist- 
ing of a free conduction band and an on-site antiferro- 
magnetic Kondo interaction which favors non-magnetic 
singlet formation. In the second order it also leads to 
the effective inter-site Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interactions which favor magnetic order. Their 
competition leads to the appearance of the quantum crit- 
ical point. In such a picture only spin degrees of freedom 
are involved in the quantum phase transition. Therefore 
the Kondo lattice model may be replaced by a simpler 
model where the itinerant hopping part is simulated by 
an inter-site interaction of the itinerant spins. Kondo 
necklace (KN) model has been originally proposed by 
Doniach^ for the one-dimensional case as a simplified ver- 
sion of the itinerant Kondo lattice (KL) model 5 . Thereby 
the kinetic energy of conduction electrons is replaced by 
an intersite exchange term. For a pure XY-type inter- 
site exchange this may be obtained by a Jordan- Wigner 
transformation. The intuitive argument is that at low 
temperatures the charge fluctuations in the Kondo lat- 
tice model are frozen out and the remaining spin fluctua- 
tion spectrum can be simulated by an antiferromagnetic 
inter-site interaction term of immobile r spins coupled 
by a Kondo interaction to the local noninteracting spins 
S. 

The quantum phase transition between the spin liq- 
uid and antiferromagnetic phases in the one dimen- 
sional Kondo necklace model has been studied by several 
methods^^. The phase diagram of the genuine Kondo 
necklace model has been studied by mean field approach 9 



in spatial dimension D=l, 2, 3. The one dimensional case 
is always in the Kondo singlet phase which is justified by 
the numerical Monte-Carlo^ and Density Matrix Renor- 
malization Group 8 (DMRG) results. While a quantum 
phase transition between the Kondo singlet phase and 
the antiferromagnetic (AF) ordered phase happens by in- 
creasing the ratio (t/J) of inter-site to the local exchange 
interaction for D=2, 3. The effect of both intersite and 
intrasite anisotropies on the quantum critical point have 
also been investigated by the same approaches. Mean 
field theory has also been applied to study the effect of 
magnetic field on the Kondo necklace modet^. However, 
for D=l the mean field approach always shows a non- 
magnetic Kondo singlet phase even in the presence of 
anisotropy in the easy axis terml*. 

In the present work we want to study the possible 
quantum phase transition in the D = 1 Kondo-necklace 
model under the assumption of the anisotropy in the 
XY-interaction (77) between the spin of itinerant elec- 
trons. However, the real space renormalization group 
(RG) studyi^ shows that if the anisotropy in the XY- 
interaction is greater than some nonzero value a phase 
transition from the Kondo singlet to the antiferromag- 
netic long range order takes place by changing the local 
exchange coupling (J/t). It means that for r] < r] c the 
system is always in Kondo singlet state but for rj > r\ c 
there is a critical exchange coupling (J c /t) which is the 
border between the Kondo singlet and AF long range or- 
dered phases. However, we will show that the quantum 
phase transition between the Kondo singlet and AF or- 
der exists for any nonzero value of r\. Our result is based 
on the spin wave approach which is justified in the ab- 
sence of local exchange for the anisotropic case. This is 
in agreement with the general statement which tells us 
about the universality of this transition. Moreover, the 
Lanczos numerical computations justify that our state- 
ment is correct. 

In Sec|TT]we will explain why the mean field approach 
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fails to produce the correct results in one spatial dimen- 
sion and easy axis anisotropy. We then consider a more 
general case of anisotropy in the Hamiltonian in Sec lIIII 
We then implement a spin wave theory to the Hamilto- 
nian with anisotropy in Sec lIVI and show that the local 
exchange will increase the quantum fluctuations which 
destroy the AF ordered state. In Sec|V]we present our nu- 
merical Lanczos results which verifies our proposed phase 
diagram. Finally, the conclusion and discussion are pre- 
sented in SeclVIl 



II. MEAN FIELD THEORY 

In this section we shortly review the mean field ap- 
proach to the Kondo-necklace model. We will show that 
this approach fails to represent the antifcrromagnetic 
phase for Ising like anisotropy. Let first consider the fol- 
lowing Hamiltonian 

Hakn = tJ2 (t*t* + rfrj + 5t*t;) + J^n ■ Si, (1) 

(ij) 

where t is the strength of exchange coupling between the 
spin (r) of itinerant electrons and J is exchange cou- 
pling between the spin of itinerant electrons and localized 
spins (S). The anisotropy parameter is represented by 
5. The bond operator representation^ is used to trans- 
form the spin operators in terms of bosonic operators*- 
(s,t x ,t y ,t z ). We start from the strong coupling limit 
where the mean field theory works well and the model 
has a Kondo singlet ground state. For the mean field 
approach^ we assume a singlet condensation (s) = s and 
a condensation for one triplet component to induce the 
antifcrromagnetic order tk, z = "vNi5k,Q+T]k,z, where r\k,z 
is the quantum fluctuation above the triplet condensation 
and Q = tt. After performing the Bogoliubov transfor- 
mation the Hamiltonian in momentum space becomes 

^AKN = + ]T[ ]T U a (k)tlj k , a + U> z {k)7)lj k , z ] 

k a=x,y 

(2) 

where Eq is the ground state energy, u) a {k) are the 
excitation energies (a = x,y,z) and $ k x{yV f)\ z (t k ,x( y ), 
fjt ) are the bosonic creation (annihilation) operators in 
the diagonal representation. (More details can be found 
in RefflQI). The excitation energies are 

w Q (fc) = V<£(*0-4/2(fc), a = x,y,z. (3) 



and 7(fc) = cos(fc). The ground state energy is 



where 



ts 2 



f x (k) = fy(k) = — 7 (fc), 

j ts 2 

d x (k) = d y (k) = u + - + —j(k), 

f z (k) = ^ 7 (fc), d z {k)=n+±+ t -£*y(k), (4) 



e£ f = e + nP 



u H ts 2 S 
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(5) 



where 



E = iv(( M -^) s - 2 - M ) + i J2 ("«(k)-d a {k) 



k,a—x,y,z 

(6) 

The chemical potential (fj,) has been added as a Lagrange 
multiplier to the mean field Hamiltonian to preserve the 
dimension of Hilbert space on each bond. The ground 
state energy should be minimized with respect to s, i 
and /x. It gives the following equations 
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J 




~ 4' 
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E(v / i+7(fc) + 2 v l 



7 (fc) 



(7) 



which should be solved self consistently. For 6 > I, we 

expect to have the long range antiferromagnetic order in 

z direction. However, the summation -i Y\ , 1 in 

N ^ k v/i+tW 

Eq. ([7]) diverges which shows that the mean field approach 
fails to work correctly in this case. 

We have even considered the extreme case of only an 
Ising term in the interaction between itinerant spins, 
namely t — > and t5 — constant. It is clear that for 
J = the ground state is long ranged antiferromagnetic 
ordered. However, the above mean field theory fails to 
show the nonzero AF ordering. We conclude that the 
mean field theory is not suitable to show the long ranged 
antiferromagnetic phase in one dimension. It should be 
related to the strong quantum fluctuations which have 
dominant effect in one dimension. Moreover, the above 
mean field theory starts from the strong coupling limit 
(t = 0) where the singlet condensation is supposed to be 
the case while for the antiferromagnetic phase we need 
to consider the weak coupling limit (J = 0). Having this 
in mind, we start from the weak coupling limit with an 
AF ordered ground state for the Ising anisotropy case. 
The effect of nonzero local exchange (J / 0) increases 
the quantum fluctuations which finally destroy the AF 
long range order. In this respect, we apply the spin wave 
theory to the anisotropic Kondo necklace chain. 



III. THE MODEL HAMILTONIAN 

The Hamiltonian defined in Eq.|T]) is always in the 
Kondo singlet phase for XY anisotropy (S < 1) while it 
is expected to have a quantum phase transition to the 
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AF ordered state at a critical exchange coupling (j|c) 
for Ising anisotropy (5 > 1). To consider a general case 
we will investigate the phase diagram of the following 
Hamiltonian in one dimension 



n = t^4 + (i - v)r!if +1 ) + ■ Si), (8) 

i=l i=l 

where rj is the anisotropy parameter which defines the 
deviation from the Ising limit (r/ = 1). The main rea- 
son to consider the anisotropy in this form is that the 
anisotropy in the itinerant part of the interaction enables 
us to investigate the effect of symmetry on the results 
of the Kondo necklace model compared with the Kondo 
model. The present model, the case of -q = has U(l) 
symmetry and for rj ^ it has Z2 symmetry. In the case 
of full anisotropy (77 = 1 , Ising limit) we expect to have 
a quantum phase transition from the antiferromagnetic 
long range order to the Kondo singlet phase. Moreover, 
in the absence of local exchange coupling (J = 0) the 
universal behavior of the XY term is the same as for all 
rj 7^ 0. It is also another motivation to see what is the 
effect of 77 anisotropy on the phase diagram of the model. 
Specially, our final results show that there exist an anti- 
ferromagnetic phase for any nonzero value of 77 when J 
is smaller than the critical exchange coupling J c (j])- The 
latter result is in contradiction with the renormalization 
group result presented in Reffl2l. 
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FIG. 1: (Color online) One dimensional Kondo-necklace lat- 
tice labeled by two sublattices. 



We then implement the Holstein-Primakov 
transformation 14 to write the spin operators in terms 
of bosonic operators. Starting from the weak coupling 
regime (J = 0) the polarized ground state is the Neel 
state which impose the following transformations for 
spins of the itinerant electrons 



r t zB = -T + aUi, 



(10) 



,+A 



( 7^ )(1 _^i)i/ 2a . 

= (V27)4a-4?) 1/2 



-A 



.-B 



(V2r)at(l- 
= (V2f)(l- 



2r 
2t 



.1/2 



a i a i -a/2 

2r ' 



In the above equations r is the magnitude of spin and the 
Boson creation (oJ) and annihilation (a) operators obey 
the following commutation relations, 



[aj,Oj] = h,i , [<4> a j] = , [oi.a.,-] = 0. 



(11) 



A similar transformations is also defined for the localized 
spin (S) by replacing r — > S and a — > e in Eqs. fTD]) 
where (e) represents the creation (annihilation) Boson 
operator for the localized spins. Generally, one should 
define two types of Boson operators (in each sublattice) 
for both itinerant and localized spins. However, due to 
the translational invariance symmetry of the model the 
Boson operators for itinerant spins in each sublattice de- 
fine the same operator as defined in Eqs. (fT0|) . The same 
story exists also for the localized ones. 

The linear spin wave theory is implemented here where 
the spin operators are described in linear form of the 
Boson operators and higher order terms have been ne- 
glected, 



IV. SPIN WAVE APPROACH 

The spin wave theory is constructed on a ground state 
which is maximally polarized. The low-lying excitations 
are expected to be created by flipping a spin and letting 
this disturbance to propagate in the crystal which is a 
spin wave. These excitations are called magnons which 
are bosonic quasi particles. 

To apply the spin wave approach the one dimensional 
Kondo-necklace lattice is labeled by two sublattices A 
and B (see FigQJ . The Hamiltonian {8| can be written in 
the following form where represents the summation 
over the nearest neighbors 

N/2 

n = *EEw^ + c 1 - 



i—l a 



N/2 



-jJ2(n A ■ Si B + H B ■ Si 



(9) 



V2t 
~2i~ 



(a, - a|) 



2 i 



(at 



(12) 



Similar expressions to Eq. (fT2"|) are applied for the other 
component of spin operators where they are replaced by 
Boson ones within the linear spin wave approximation 
which finally leads to the following form for the whole 
Hamiltonian 



n 



N{tr 2 - 

N/2 



N/2 

Jr^)+tr^^(2 



i=l 



2 

N/2 



N/2 



reUi 



2jVrSj2(a 1 e l + a\el). (13) 



To proceed further and diagonalize the Hamiltonian we 
first transform the operators to their Fourier counterparts 
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by the following relations, 



Cfc 



(14) 



where the ground state energy is 



E (J) = N[--(l + J) + VT+2J], (20) 



fc j 

(15) 

In the momentum space representation and for r = s = 
1/2 the Hamiltonian (Eq. (fT3|) ) is finally written in the 
following form 

H = ~(l + J) + ^-^ Yl 7(A)[4ct fc + ^fc] 



and the quasiparticle excitations are 



fc 



k 

+4 h -k + h\h k + hl k h- k ) 

+ £[(.7/2 + 1) - ii_^ 7 (fc)] [c t Cfc + cL fc c_ fe ]. (16) 
fc 

The Hamiltonian can be diagonalized by the para-unitary 
transformation— for the general case of r\ ^ 1. However, 
for the case of r\ = 1, one can use the Bogoliubov transfor- 
mation which is a special case of the general para-unitary 
transformation. 



The Hamiltonian of the general case ij ^ 1 can be 
diagonalized by the para-unitary transformation^ for the 
bosonic operators. In this formalism the Hamiltonian is 
written as a hermitian 2m-square matrix D. The para- 
unitary transformation (r) can be constructed if and only 
if D is positive definite. This procedure is presented in 
appendix. 

However, for the special case r\ = 1 the para-unitary 
transformation is the usual Boguliubov transformation 16 . 
In this case, the bosonic operators is transformed to a 
new set of bosonic operators (a, (3) by the following rela- 
tions, 



c k = u k Oik + Vkfl'lk , hk = u k (3 k + v k o?_ k 



(17) 



where u k = cosh(0/c) and v k = sinh((9fc) and 9 k is the free 
parameter of transformation. The parameter 8 k should 
satisfy the following relation to have a diagonal Hamilto- 
nian 



tanh(26» fe ) 



j + r 



(18) 



which is actually independent of k. The diagonalized 
Hamiltonian in terms of the new sets of operators is 



W a {k) = -1 + V1 + 2J , uSp(k) = 1+V1 + 2J. (21) 

The sublattice magnetization itia — jjJ2i( T ~ a \ a i) 
defines the order parameter for the AF long range ordered 
phase which can be calculated in the diagonal bases of 
the Hamiltonian. For simplicity we put t — 1 as the scale 
of energy. In the weak coupling limit J = the Hamil- 
tonian (Eq.©) has long range AF order. The quantum 
fluctuations which are created by nonzero J decrease the 
amount of AF ordering and finally destroy the ordered 
phase at a critical value, J c , where the sublattice magne- 
tization becomes zero 



1 / t v 



= 1 - 



■7 + 1 
2V1 + 2J 



= 



J r ~ 6.46. 



(22) 

It is the first estimate of the critical value of the exchange 
coupling. However, the condition which we have stated is 
the extreme value where the quantum fluctuations sup- 
press the sublattice magnetization totally. At this ex- 
treme condition the accuracy of the linear spin wave re- 
sults is not justified, since the amplitude of the quantum 
fluctuations is not small. However, this approach is valid 
for small amplitude of fluctuations. Moreover, it shows 
that the amplitude of fluctuations become large by in- 
creasing the value of exchange coupling (J). Thus, the 
qualitative picture of the spin wave approximation is cor- 
rect. 

To justify that the spin wave theory gives correct re- 
sults, we have plotted the ground state energy per site 
versus the exchange coupling (J) in Fig(2j We have com- 
pared the results obtained in Eq. (|2"0|) with the numerical 
Lanczos ones which will be discussed in the next section. 
It shows very good agreement of the spin wave with Lanc- 
zos results which justifies that the spin wave ground state 
presented here can be a good candidate to represent the 
many body ground state in this case. 



B. 



We have implemented the para-unitary transformation 
(see appendix) to diagonalize the Hamiltonian (Eq.(Ql)])). 
The diagonalized Hamiltonian in terms of new sets of 
bosonic operators (<j), ip) is 



H = E o {J)+^(W a {k)ala k +uJ (k)/3ll3k) (19) H = E Q ( V ,J) + ^[fi Q (fc)0^ fc + n fJ (k)^ k ]. (23) 

fc 
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j 



FIG. 2: (Color online) Ground state energy per site 
(Eo(J)/N) versus J for 77 = 1 obtained by numerical Lanczos 
method (extrapolated to the N — > 00 with 4 digits accuracy) 
and linear spin wave theory which show very good agreement. 



where Ea(r),J) = (1/N) J2k e o(Vi J> k) 1S the ground 
state energy and eo(r?, J; k) has the following expression 

e ( v ,J;k) = ^= {VW~~Q + VvTQ) , 

V = 4(1 + J) - (2 + J) 7 (fc) , Q = V^ + T 2 , 
5 = (2 + J) 7 (fc)»? - 8J 2 (2 + 7 (A)(i, - 1)), 
T = 4 + J(4 + 7 (*)(t7 - 1)) + 2 7 (fc)(77 - 1). (24) 

The frequencies r2 a (fc) and Qp(k) are 



n a (fc) = -7=V^(i + J) - (2 + j)j{k) - a 



Qp(k) = — y/4(l + J)- (2 + jyy(k)+ A, (25) 
where 

A = VS + T 2 . (26) 

A remark is in order here, the results in this subsection 
will not reproduce the results of 77 = 1 case simply by 
putting r\ equal to one. The reason is related to a diver- 
gent factor of (1 — 77) _1 which appears in the intermediate 
part of calculations. Such terms will not appear if we put 
X] = 1 from the beginning and the calculations will be- 
come much simpler as presented in the last subsection. 
In the diagonal bases of the Hamiltonian we can calcu- 
late the sublattice magnetization as defined in Eq. (fTQ|) . 
As we have discussed in the previous subsection we can 
determine the critical border between the AF order and 
the Kondo singlet phase by looking for the position where 
the sublattice magnetization becomes zero. 

Our results which are summerized in TableHJ show 
that the onset of nonzero XY-anisotropy (77 ^ 0) develops 
the long range AF order for small local exchange coupling 
(J < J c ). This is also justified by the numerical Lanczos 
results where the static structure factor at momentum tt 



diverges for any nonzero 77 and J < J c . However, the 
numerical values for J c which are obtained by the spin 
wave approach is far from the numerical Lanczos ones. 
This will be discussed in the Sec lVIl 

V. NUMERICAL LANCZOS RESULTS 

The numerical Lanczos method has been applied to 
the one dimensional anisotropic Kondo-necklace model 
defined in Eq.JS]). The anisotropic Hamiltonian does not 
commute with the total z-component spin. Thus, one 
should consider the full Hilbert space for doing numer- 
ics. In this respect, we have performed the numerical 
computations with TV = 12, 16, 20, 24 spins. The com- 
putations have been done for different values of the ex- 
change coupling (J) and the anisotropy parameter (77) to 
trace the quantum phase transition between the Kondo 
singlet phase and the antiferromagnetic one. The Kondo 
singlet ground state is not ordered while the antiferro- 
magnetic phase has long range order. The long range 
order imposes that the structure factor diverges in the 
thermodynamic limit (N — » 00) for the antiferromag- 
netic wave vector k — n. The structure factor of the 
z-component spin at the wave vector k is defined by 




k 



FIG. 3: (Color online) The z-component structure fac- 
tor versus k for r\ = 0.6, N = 24 and different J = 
0.0,0.1,0.2, ...,1.0. 

We have plotted the structure factor (G(k)) versus the 
wave vector (fe) for 77 = 0.6, N — 24 and different ex- 
change couplings (J = 0.0, 0.1, 0.2, . . . , 1.0) in Fig(3] The 
structure factor gets a strong peak at k = tt for some 
values of J. This peak shows the presence of antiferro- 
magnetic order. We will show in the next plots how the 
height of this peak evolve as the size of system changes. 

The z-component spin of structure factor at the an- 
tiferromagnetic wave vector, k — 7T, has been plotted 
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FIG. 4: (Color online) The z-component structure factor for 
k = tt versus the inverse of of system size (1/N) for 77 — 0.0 
and different J = 0.0, 0.1, 0.2, . . . , 1.0. The structure factor 
diverges for J = like y/N which represent the massless be- 
havior of XX spin 1/2 chain. The other sets of data show no 
divergence or weaker than y/N which justifies no long range 
antiferromagnetic order. 



in Figs|4][8] for different anisotropy (77) versus the in- 
verse of system size (N). The system sizes are N — 
12, 16, 20, 24. Each figure contains different plots for 
J = 0.0, 0.1, 0.2,..., 1.0. 

We have plotted the case of genuine Kondo-necklace 
model, 77 = in FigJU For J = 0, the properties of the 
XX spin 1/2 chain is reproduced which shows that the 
structure factor diverges^ proportional to VN. How- 
ever, it does not show the antiferromagnetic order of this 
model. It only represent the critical properties of the 
XX spin 1/2 chain which has algebraic decay of correla- 
tion functions. Turning on the local exchange to nonzero 
values J /: we observe weaker growth of the structure 
factor as the system size in increased. It is concluded that 
the one dimensional genuine Kondo-necklace model does 
not have an antiferromagnetic phase and is always in the 
Kondo singlet state. Moreover, we expect the structure 
factor diverges stronger than \/N for the antiferromag- 
netic ordered phase. 

The presence of anisotropy (77 7- 0) breaks the U(l) 
symmetry and reduces it to Z2 symmetry. This is the 
onset of antiferromagnetic ordering formation which is 
the result of symmetry breaking. The structure factor at 
k = tt versus the inverse of length (1/N) for rj = 0.1 has 
been plotted in FigJS] It is clear that G(ir) diverges faster 
than \f~N for J = which verifies the formation of an- 
tiferromagnetic ordering in the ground state. A similar 
situation exists also for J — 0.1,0.2,0.3,0.4. However, 
for J > 0.5 the divergence of G(ir) is not apparent as 
N — > 00. We then conclude that for J < J c (rj = 0.1) ~ 
0.4 the model has antiferromagnetic ordering while for 
larger values of J it is in the Kondo singlet phase. Thus, 
a quantum phase transition occurs at J c from the antifer- 
romagnetic order to the disordered Kondo singlet state. 

We have also plotted the z-component structure factor 



FIG. 5: (Color online) The z-component structure factor ver- 
sus 1/N (inverse of size of system) at k = tt for 77 = 0.1 and 
different J = 0.0, 0.1, 0.2, . . . , 1.0. The divergence of G(k = tt) 
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FIG. 6: (Color online) G(k — n) (the z-component struc- 
ture factor) versus 1/N for 77 = 0.2 and different J = 
0.0,0.1,0.2, 1.0. It is clear that G(k = tt) diverges for 
J < 0.4 as N goes to infinity. 



at antiferromagnetic wave length (G(k — tt)) versus the 
inverse of system size (1/N) for 77 = 0.2 in FigJBl 77 = 0.6 
in Fig[7| and 77 = 1.0 in FigfS] We have observed that 
G(k = tt) diverges in the thermodynamic limit (N — * 00) 
for J < J c (r/) which verifies the existence of antiferromag- 
netic ordering. Our results for the critical exchange have 
been summerized in table |U The accuracy of our results 
for the critical point (J c ) is about ±0.05. Because, the 
size of system is limited to JV = 24. To distinguish be- 
tween the ordered and disordered phase more accurately 
it is necessary to consider larger system sizes. 



VI. CONCLUSION AND DISCUSSION 

We have studied the phase diagram of the one dimen- 
sional anisotropic Kondo necklace model. We have fo- 



TABLE I: The critical exchange value (J c ) for different values 
of anisotropy (77). The Lanczos results are read from the data 
of G(n) versus 1/N which has an accuracy about ±0.05. 
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FIG. 7: (Color online) The z-component structure factor 
versus 1/N for k = n and r\ — 0.6 and different J = 
0.0, 0.1, 0.2, . . . , 1.0. The structure factor diverges for J < 0.4 



cused our attention on the case of Ising anisotropy (easy 
axis of the spin of itinerant electrons, S > 1). Let us sup- 
pose that the interaction between the spin of itinerant 
electrons is only composed of rfrj type. It is clear that 
in the absence of local exchange interaction (J = 0) the 
model has a classical antiferromagnetic ordered ground 
state with nonzero staggered magnetization. It is ex- 
pected that a nonzero local exchange (J ^ 0) adds quan- 
tum fluctuations which destroy the antiferromagnetic or- 
der at a critical exchange (J c ). However, the mean field 
approach^ fails to produce this picture correctly. In 
the mean field approach which is based on the strong 
coupling limit (J/t — ► oo) we have assumed a nonzero 
condensation for local singlet formation plus a nonzero 
triplet occupation (f) which induces the antiferromag- 
netic order. The solution of the self consistent mean 
field equations always give no magnetic order, i = 0. 
A divergent integral appears if we assume t ^ 0. How- 
ever, the mean filed approach works fairly well in two and 
three dimensional models in addition to the case of XY 
anisotropy (0 < S < 1) of the one dimensional system^. 

To overcome the problem which appears in the mean 
field approach for Ising anisotropy of the one dimensional 
model we decided to start from the weak coupling limit 
(J — * 0). The linear spin wave theory has been imple- 
mented which is based on a classical antiferromagnetic 
ground state. To study the phase diagram we have con- 
sidered a general case of Hamiltonian which has been 
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FIG. 8: (Color online) The z-component structure factor 
versus 1/N at k = n for r\ — 1.0 and different J = 
0.0, 0.1, 0.2, . . . , 1.0. The divergence of G(k = tt) for N -> oo 
is clear for J < 0.5. 

given in Eq.([8|). Our results show that the quantum fluc- 
tuations grow with increasing J such that the sublat- 
tice magnetization becomes zero at a critical exchange 
coupling (J c )- Moreover, the antiferromagnetic phase 
exists for any nonzero value of anisotropy (rj ^ 0) and 
J < J c . It is in contradiction with the results presented 
m 

RefG3 which gives the antiferromagnetic phase only 
for -q > 7] c ~ 0.58. In other words our results state that 
rjc = 0. We have plotted the phase diagram of the one di- 
mensional anisotropic Kondo necklace model in the J — rj 
plane in Fig|9] It is obvious that for any nonzero rj there 
exists an antiferromagnetic phase. Moreover, the critical 
exchange coupling (J c ) depends strongly on rj for small 
anisotropy, < rj < 0.1 while its dependence is weak for 
rj > 0.1. 

Our result is in agreement with the symmetry argu- 
ments which is based on a qualitative picture. Let us 
first suppose that J = 0. The remaining Hamiltonian is 
an anisotropic XY spin 1/2 chain. The quantum renor- 
malization group^i verifies that the universality class of 
the nonzero anisotropy (77 7^ 0) is the same as Ising case 
(77 = 1). Thus, if the antiferromagnetic phase exists for 
the Ising case it should also appear for any nonzero value 
of < rj < 1. Moreover, the Hamiltonian has U(l) sym- 
metry at rj = while the symmetry is broken to the lower 
Z2 for rj 7^ 0. It is thus expected that the quantum phase 
transition happens at the symmetry breaking point. We 
expect that the addition of local exchange (J 7^ 0) be- 
comes irrelevant at the fixed point J = 0, which is in 
agreement with our numerical results. However, this 
needs to be proved by more sophisticated methods like 
weak-coupling renormalization group which is out of the 
scope of this article. 

To verify that our result gives the correct phase 
diagram we have implemented the numerical Lanczos 
method. Our results show that the z-component struc- 
ture factor diverges at the antiferromagnetic wave vector 
for any nonzero rj and J < J c . However, the value of 
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FIG. 9: Phase diagram of the one dimensional anisotropic 
Kondo necklace model where a quantum phase transition 
(dashed line) separates the ordered antiferromagnetic phase 
from the Kondo singlet one. 



the quantum critical point (J c ) can not be determined 
with very high accuracy because of the finite size effect. 
In this study we have been limited to N = 24 because 
the full Hilbert space should be considered in computa- 
tion and the Lanczos is basically limited to such values. 
However, it is proposed that a density matrix renormal- 
ization group can get more accurate values. Moreover, we 
have understood that the critical exchange ( J c ) which is 
given by spin wave theory is far from the real one which 
is given by numerical Lanczos method. The reason is re- 
lated to the strong quantum fluctuations which exist close 
to critical point and ruin the accuracy of linear spin wave 
results. 
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APPENDIX A: THE PARA-UNITARY 
TRANSFORMATION 

This section is based on the approach proposed in 
Reffl5l to diagonalize a bosonic Hamiltonian. The Hamil- 
tonian in Eq. (fT5|) can be written in the following matrix 
form 



where is a row vector 



rt = (4 h l c -k h- k ), 



(Al) 



(A2) 



and D is a square matrix 

7 + 4-(l-r7)7(fc) 



D = 



V 





(l-r?b(fc) 
4 
J 
2 



(l-r7)7(*0 
4 

L 1 

J 2J+4-(l-77) 7 (fc) 
2 4 








i) 



(A3) 



The Hamiltonian will be diagonalized by the following 
transformation (J) to the new set of boson operators, 



where 



A = JT, 



At = (4 <]>_ k V- fc ), 



(A4) 



(A5) 



and the diagonal representation is given in Eq . (123)) . The 
new set of operators should satisfy the bosonic commu- 
tation relation, namely 

[Aj, A]] = Si j, (A6) 
where 5 is the para Kronecker symbol defined by 



Si,i = 1 
Si.i = -1 



J i,3 







1 < i < 2, 
3 < i < 4, 
i 3- 



(A7) 



In other words the transformation matrix (^7) should 
obey the following equations, 



jiji = i or j^ij = i, 



(A8) 



where I is the para unit matrix, / = diaqjl ,1,-1,-1). A 
theorem which has been proved in Ref J 1 51 states that the 
matrix D can be para unitary diagonalized into a matrix 
with all diagonal elements positive if and only if D is 
positive definite. According to Eq. (|A3p . D is positive 
definite except at the point (rj = Q, J = 0). To find the 
para unitary transformation which satisfies Eq. (|A8|) the 
Hamiltonian and the transformation will be written in 
the following matrix form 



D = 



A B 
B A 



and 



J 



-ut rjt 



(A9) 



(A10) 



where A,B,U and V are 2x2 matrices. 

The 2x2 matrix K is defined by K^K = A - B. 
The eigenvalues of the Hermitian matrix K(A + B)K^ 
are identified by det[K(A + B)K^ - A? 7] = 0, i = 1, 2 
where / is the 2x2 unit matrix and the corresponding 
eigenvectors (xi) are normalized such that xlxi 
We define fi — K^Xi an d obtain pi such that 



A, ' 



(A + B)fi = XiPi. 



(All) 
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The transformation matrices U and V are given by their REFERENCES 
column vectors with the following relations, 

V>i = \{fi +Pi), 

Vi = l(fi-Pi). (A12) 



* langari@sharif.edu; http://spin.cscm.ir 

1 M. A. Continentino, Brazilian Journal of Physics 35, 197 
(2005). 

2 M. Vojta, Rep. Prog. Phys. 66, 2069 (2003). 

3 P. Fulde, P. Thalmeier and G. Zwicknagl Solid State 
Physics Vol. 60 (Academic Press, 2006). 

4 S. Doniach, Physica 91B, 231 (1977). 

5 H. Tsunetsugu, M. Sigrist and K. Ueda, Rev. Mod. Phys. 
69, 809 (1997). 

6 R. Jullien, J. N. Fields and S. Doniach, Phys. Rev. B 16, 
4889 (1977). 

7 R. T. Scalettar, D. J. Scalapino, and R. L. Sugar, Phys. 
Rev. B 31, 7316 (1985) and references therein. 

8 S. Moukouri, L. G. Caron, C. Bourbonnais, and L. Hubert, 
Phys. Rev. B 51, 15920 (1995). 

9 G. M. Zhang, Q. Gu and L. Yu, Phys. Rev. B 62, 69 (2000). 



A. Langari and P. Thalmeier, Phys. Rev. B 74, 24431 

(2006) . 

1 P. Thalmeier and A. Langari, Phys. Rev. B 75, 174426 

(2007) . 

2 A. Saguia, T. G. Rappoport, B. Boechat and M. A. Con- 
tinentino, Physica A 344, 644 (2004). 

3 S. Sachdev and R. N. Bhatt, Phys. Rev. B 41, 9323 (1990). 

4 T. Holstein and H. Primakoff, Phys. Rev. 58, 1908 (1940). 

5 J. H. P. Colpa, Physica A 93, 327 (1978). 

6 N. Bogoliubov, J. Phys. (USSR) 11, 23 (1947), also 
reprinted in D. Pines, The Many-Body Problem (W. A. 
Benjamin, Inc., New York, 1961) p. 292. 

7 A. Langari, Phys. Rev. B 69, 100402(R) (2004). 

8 M. Karbach, K.-H. Mutter and M. Schmidt, Phys. Rev. B 
50, 9281 (1994). 



